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Abstract 

We present a consistent analysis of final state interactions in K — > 3tt de- 
cays in the framework of Chiral Perturbation Theory. The result is that the 
kinematical dependence of the rescattering phases cannot be neglected. The 
possibility of extracting the phase shifts from future K$ — Kl interference 
experiments is also analyzed. 
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I. INTRODUCTION 



Unitarity requires that final state strong interactions should be taken into account in 
kaon decay amplitudes. On the other hand, such effects are interesting by themselves, as 
they reflect the properties of low-energy meson interactions. In this regard, they could 
represent a meaningful test of the theoretical approach based on effective chiral Lagrangians 
and the related amplitude expansions in powers of meson momenta (Chiral Perturbation 
Theory), which should incorporate general features of long-distance QCD [1|,0]. 

The role of n — n phase shifts in K — > 2n decays has been extensively discussed [|J. 
These transitions involve the 7r — 7r phase shifts at y/s = tuk, and at this energy scale the 
correction to the leading order prediction in Chiral Perturbation Theory (ChPT) are found 
to be large in each isospin channel Hj. 

In K — > 3tt, final state strong interactions act at substantially lower energy compared to 
K —>■ 2tt, and consequently in current fits to experimental data they have been assumed to 
be negligible |5|U. As a matter of fact, these effects could be experimentally accessible in 
the near future. 

In the approximation of considering only two-body strong interactions, thus neglecting 
the irreducible 3n rescattering diagrams which should be suppressed by phase space, in 
principle K — > 3ir final state interactions can provide a complementary way to study n — n 
phase shifts near threshold. Indeed, this is the region where these phase shifts can be most 
reliably predicted in the framework of ChPT. Therefore, such an analysis should allow an 
independent, and significant test of this theoretical approach, to be combined with, e.g., 
those from decays @, also relevant to the n — ir low-energy region. 

As another important point of interest, we recall that the knowledge of rescattering is 



crucial in order to estimate direct CP- violating asymmetries in K — > 3tc [p|HTl||. Manifesta- 
tions of such asymmetries would allow to determine the existence of direct CP-violation in 
a channel alternative to K — > 2tt, and to improve our knowledge of this phenomenon, which 
is predicted by the Standard Model but not clearly established yet. 



As emphasized in [12], a convenient way to study final state interactions in K — > 3n is to 



consider Kl — K$ interference in vacuum as a function of time at "interferometry" machines 



such as LEAR [13[ and the 0-factory DA$NE p^Jlq] . The typical interference term has the 
form 

[(&k\K s )*{&k\K l ) exp (iAmt)] exp (- S * L t), (1) 

where Am = rriL—ms is the K^—Ks mass difference. Indeed, studies of the time dependence 
of Eq. (HP should lead to a determination of both the real and the (expected small) imaginary 
parts of the amplitudes. The advantage is that the latter ones appear linearly in ([!]), whereas 
they appear quadratically and consequently are much less accessible in width measurements. 
To obtain an order of magnitude estimate of the effect of the interference term, constant 
phase shifts (independent of energy) were assumed in Ref. fl2f . 

Actually, in the case of K — * 3tc it is not quite appropriate to use the notion of constant 
phase shifts, because: i) there are two independent 1=1 final states which can be connected 
by strong interactions, so that one should introduce a 2 x 2 mixing (or rescattering) matrix; 
and ii) in general the rescattering matrix elements are functions of pions momenta. 
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Phenomenologically, K — > 3tt transition amplitudes for the various modes are expanded 
in powers of the kinematical variables with constant coefficients, so that the momentum 
dependence of rescattering should be taken into account for a consistent expansion. This 
is particularly desirable also in connection with momentum expansions predicted by ChPT. 
Previous theoretical estimates of momentum-dependent K — ► 3n strong rescattering were 
performed in the non-relativistic approximation in Ref. |l6j and in leading order ChPT for 
charged kaons in Refs. |8| and || . In the framework of a complete analysis of all the ChPT 
quartic effects in K —>■ 3tt, the imaginary parts were expanded in the kinematical variables 
and calculated numerically in Ref. |J. 

In this paper, we extend the calculation of Refs. || and || to obtain the K° — > 3n 
rescattering matrix in ChPT. To this purpose, we shall first review the general symmetry 
and unitarity constraints on K — > 3ir amplitudes, then we will construct a convenient, model 
independent parametrization of the strong rescattering matrix which is unique for all decay 
modes, and finally we will use ChPT to estimate it. The results can be applied to make 
more reliable predictions for the K ^ 3ir time correlations. In addition, our analysis allows 
to clarify some delicate questions regarding direct CP- violation in K — > (37:^. 

Specifically, the plan of the paper is as follows: in sect. Jljwe set the formalism to expand 
K — > 3tt amplitudes; in sect. [ITJ we define the rescattering matrix and evaluate it in lowest 
order ChPT; in sect. ^ we discuss some consequences for the Dalitz plot analysis and 
CP- violation; in sect. [V] we present the expectations for the interference term; finally, sect. 
\/T\ contains some concluding remarks. 



II. K -> 3vr FORMALISM 
In the limit of CP conservation, there are five different channels for K — > 3tt decays: 



K ± - 


-> 7r ± 7T ± 7T =F 


(1 = 


1,2 


K ± - 


-> 7r ± 7r°7r° 


(1 = 


1,2 


K L - 


-> 7T ± 7r =F 7r° 


(/ = 


1) 


K L - 


-> ttVtt 


(/ = 


1) 


K s - 




(/ = 


2). 



(2) 



Here in parentheses we indicate the isospin values relevant to the final (3tt) states, assuming 
only AI = 1/2, 3/2 transitions. In principle, the Ks decay to the 1 = state is not forbidden, 
but due to Bose symmetry it is strongly suppressed by a high angular momentum barrier 



17| and we neglect it. The first four modes are dominated by AI =1/2 transitions, while 



the last one is a pure AI = 3/2 transition and only recently has become accessible through 



time-dependent interference experiments [18 



For K(p) — ► 7ri(j9i)7r2(p2)7T3(P3) decays we introduce the familiar kinematical invariants 

Si = (pK-Pi) 2 and s Q = )-^Si = )-m 2 K + ml, (3) 

6 i 6 

where the index i = 3 refers to the "odd" charge pion. Neglecting isospin breaking effects, 
following e.g. Refs. |17|,[l9|, we can decompose the decay amplitudes in the general form 
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-4++- = 2A c (s 1 , s 2 , S3) + B c (si, s 2 , s 3 ) + B 2 (si, s 2 , s 3 ) 

A + oo = A c (si, s 2 , s 3 ) - S c (si, s 2 , s 3 ) + £ 2 (si, s 2 , s 3 ) 

A+_ = ^n(si,s 2 ,s 3 ) - B n (sx,s 2 , s 3 ) (4) 

A^ 00 = 3j4 n (si,s 2 ,s 3 ) 

^+_o = B 2 (s u s 2 ,s 3 ). 

Here, reflecting Bose symmetry and the assumed CP conservation, all amplitudes Aj and 
Bj (j = c, n, 2) are symmetric under exchange (1 <-> 2). Furthermore, the amplitudes A,- 
are completely symmetric for any permutation of the indices 1,2 and 3. Conversely, the 
amplitudes Bj do not have this symmetry, and under permutations of indices only obey the 
relation 



Bj(si, s 2 , s 3 ) + Bj(s 3 , s 2 , si) + Bj(si, s 3 , s 2 ) = 0. (5) 

Finally, the amplitude B 2 is antisymmetric for the exchange (1 2). It is not independent 
from the other ones, and can be expressed in terms of B 2 as 

2 

B 2 (s 1 , s 2 , s 3 ) = - [£2(53, s 2 , si) - B 2 (si, s 3 , s 2 )] . (6) 

Concerning isotopic spin, the amplitudes Aj and (B c ,B n ) correspond to AI = 1/2 and 
AI = 3/2 transitions to the 1 = 1 final three-pion state, while B 2 is associated to the 
AI = 3/2 transition to / = 2. 

From the decomposition above we note that there are two amplitudes leading to / = 1 
final states, which differ for the pion exchange symmetry properties, namely the A's are fully 
symmetric whereas the B's have mixed symmetry. Accordingly, it is convenient to introduce 
the two matrices 

which in the 1 = 1 sector transform the symmetric and non-symmetric amplitudes into the 
physical ones for charged and neutral kaons, respectively. Thus: 

\ A% )- Tc { BM) ) { ^00 ) ~ Tn \ B n (s t ) ) ■ (8) 
Defining the two dimensionless Dalitz plot variables 

Y='^ and X = ^, (9) 
mt mi 

TV TV 

and taking into account the symmetry properties of A's and S's, we can expand the five 
independent amplitudes in Eq. ([|) in powers of X and Y up to quadratic terms: 

Aj = aj + Cj (Y 2 + X 2 /3) 

B^bjY + d^Y'-X 2 /!). (10) 
Substituting Eq. (|i~Cf) in Eq. (f|), we obtain 
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= 2a c + (b c + b 2 )Y + 2c c (Y 2 + X 2 /3) + (d c + d 2 )(Y 2 - X 2 /3) 
A +00 = a c - (b c - b 2 )Y + c c (Y 2 + X 2 /3) - (d c - d 2 )(Y 2 - X 2 /3) 
A L + _ = a n - b n Y + c n (Y 2 + X 2 /3) - d n (Y 2 - X 2 /3) (11) 
^ooo = 3a„ + 3c n (Y 2 + X 2 /3) 
A s + _ = p 2 X - ±d 2 XY. 

This decomposition can be easily related to the one introduced in Refs. [§,|J. 



III. 3vr FINAL STATE INTERACTION 

Since strong interactions are expected to mix the two 1 = 1 final states, we must introduce 
a strong interaction rescattering matrix which mixes the corresponding decay amplitudes. 
Projecting the final state (37r)/ =1 by means of the matrices T" 1 and T~ x in the symmetric- 
nonsymmetric basis, we can define the scattering matrix R, common to charged and neutral 
channels, as follows: 



4(1) 
^+00 



R 




t c rt: 



A 
A 



(i> 



(i) 
-oo 



(12) 



^ooo 



T„R 



R 



A lt 



T n RT n 



^ooo 



(13) 



Here the subscript R means that in the decay amplitude rescattering has been included. 
The matrix R defined above has diagonal elements which preserve the symmetry properties 
under pion exchanges, as well as off-diagonal elements which connect symmetric amplitudes 
to non-symmetric ones and viceversa. 

The unitarity conditions are obtained by imposing conservation of probability, namely: 



1(41 



)R 



+00)R\ 



o)r\ 2 + |(^ooo)j?l 



\A 



(i) 



\A 



(1) 12 
+00 1 



+ \A 



000 I 



(14) 
(15) 



where <i$ represents the phase space element. 

We now perform the calculation of R using ChPT. At the lowest order p 2 , there are no 
quadratic terms and the coefficients a,j, bj in fllTf ) are real if CP is conserved. At order p 4 
loops and counterterms will appear, generating real parts with higher powers in X and Y 
and also imaginary parts proportional to the 0(p 2 ) constants aj and bj . These imaginary 
parts define the rescattering matrix relevant to the constant and linear terms. In principle, 
imaginary parts of quadratic terms can occur similarly, but since these ones appear at the 
higher order p 6 we neglect them, as it is also justified by the smallness of the experimental 
values of quadratic slopes. In this approximation, we replace A c>n and B c n in ( |T2"D and ([13]) 
by a c ,n and b c ^ n Y respectively. The unitarity conditions resulting from (|i~4|) and (|i~5|) are 
equivalent, and take the form 
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\Rn 



+ r\R2l\ 

5 



\R 



22 



+ 5 lft>l 



d<$> 

d$Y 2 



(16) 

Y' 1 / r/<I> V ' J (17) 
/" d$ [5i2ni2J 2 + 2R 21 R* 22 ] Y = 0. (18) 
As anticipated, at 0(p 2 ) there are only tree diagrams that can be easily computed with 



the leading order chiral weak Lagrangian 



and there is no final state interaction so 



that R = I (trivial case). At order p A loop diagrams (Fig.l) generate imaginary parts, 
corresponding to on-shell propagators in internal lines, so that we can write 



R = I + i 

Using the strong 0(p 2 ) chiral Lagrangian 

F 2 



a{Si) 
a'(si) 



(19) 



(2) 



-tr 



d^d^ + M(S + S f 



(20) 



where £ = exp(z0/v / 2i 71 7 r), 4> is the octect matrix of the pseudoscalar fields, F w is the pion 
decay constant (F n ~ 93 MeV) and M = diag(m 2 , m 2 , 2m 2 K — m 2 ), we find 



a{Si) 
a'(si) 



WirF* & 3 
32ttF2 ^ 3 L 



Vi[Si - ml) -v 3 (s 3 - ml) 



32ttF2 
1 



32ttF2 ^ 3 



-J2 v i( s * ~ 4ml) + J2 m l 
6 i=i »=i 

^ 2 (si-m^)(s -Si) 



2 17 3 (S 3 - S ) - Ui(Sj - S ) 



s 3 — s 



(21) 
(22) 
(23) 
(24) 



where Vi are the "velocities" : Vi 
(|Hf) is nontrivial and implies 



s 3 ~ s o 

1 — Ami/ Si) 1 ^ 2 . At this order, only the unitarity condition 



(25) 



This condition is exactly verified by the functions in Eqs. ( f22]) and Eq. (|2~4]) , as expected 
since ChPT is an effective field theory where unitarity is perturbatively satisfied. 

We remark that the rescattering matrix R could have been directly evaluated by just 
integrating the n — n scattering amplitude over the phase space of intermediate particles. 
Actually, once the matrix R has been defined, one could improve the lowest order Eqs. 
(Enp- fl24D by including all higher orders in strong interactions, or even by replacing them by 
any other available phenomenological information on it — tt scattering. Analogously, for the 
weak amplitudes a c , b c and a n , b n one could use either the ChPT predictions or the available 
experimental determinations. 
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For the decay into the 1 = 2 final states there is only one amplitude, with definite 
symmetry under pion exchange which must be preserved by strong interactions. Thus, we 
can write 



{B 2 ) R = b 2 Y{l + i5{ Si }) 
(B 2 ) R =~b 2 X(l + i5( Si )), 



(26) 



where two functions 5 and 5 are again not independent because (B 2 ) R and (B 2 ) R must satisfy 
Eq. ffl). At the lowest non trivial order in ChPT we find 



S(si 



32-kFI 



1 3 

° i=i 



'si - Aml s 



Vifa - s )(2si - 5ml) - v 3 (s 3 - s )(2s 3 - 5ml) 



s 3 — s 



(27) 



Regarding three-body rescattering, which will appear at two loops, we would expect its 
contribution to the functions a, a', j3, [3' and 5 to be rather small, as being suppressed by 
phase space, assuming that the three-body coupling is not anomalously large. This is indeed 
the case for the leading order Lagrangian fl2"0"|). Regarding higher orders in ChPT, 0(p A ) 
contributions to R might be relevant, similar to the case of ir — n phase shifts where the 
scattering lengths turn out to be affected at the 20% — 30% level j|. We can take these 
figures as an indication for the accuracy of the subsequent applications of Eqs. (^l|)-(|27|). 



IV. CONSEQUENCES FOR DALITZ PLOT ANALYSIS AND CP- VIOLATION 

As expected from the smallness of the available phase space, the functions a, a', j3, (3' 
and S are smaller than unity over the whole Dalitz plot. Indeed, by expanding in powers of 
X and Y up to quadratic terms, we obtain 

a(X,Y) ^ a + a 1 (Y 2 + X 2 /3) [ a Q ~ 0.13 a x ~ -2.9 x 10~ 3 

a'(X, Y) ^ a' Y + a[(Y 2 - X 2 /3) [ a' ~ -0.12 a[ ~ 3.4 x 10" 3 



P(X,Y) = (3 Q + (3 X (Y 2 - X 2 /3)/Y [ fa ~ 0.047 A ~ 4.7 x lO" 3 

(3'(X, Y) = {3' (Y 2 + X 2 /3)/Y [ (3' = a' /5 

5(X, Y) = 5 + 5i(F 2 - X 2 /3)/Y [ 5 = -f3 5 1 ~ -0.020 



(28) 



Using Eqs. (|i~2"D , (|13D , ( p6|) and (|28|) we can expand both real and imaginary parts of all 
K — ► 3tt amplitudes up to linear terms. 

As a first application of the formalism we can discuss the role of rescattering in CP- 
odd charge asymmetries in K ± — > 3n. Contrary to K — ► 2tt, where direct CP-violation is 
suppressed by the smallness of the AI = 3/2 amplitude, in K — > 37r an observable effect 
can potentially arise also from the interference of the two AI = 1/2 amplitudes. For 
a nonvanishing effect it is crucial that the relevant amplitudes have different electroweak 
phases (which can be the case only at order p A in the framework of ChPT) as well as 
different rescattering phases. 
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Let us consider, for example, the amplitudes for K + — > (tt + 7t + 7t )i=\. From the preceding 
relations we easily obtain 

$le(A { + } + _) = 2a c + b c Y 

Qm(A ( ^ + _) = 2a c a + a c a' Y + b c (3 Y = 2a c a + b c Y (fi + y®'o^) 



(29) 
(30) 



In Eq. (|3CiD the contribution of a' is multiplied by the sizable factor \a c /b c \ ~ 3.5 — 4.0, 
and dominates over the one of (3q by almost one order of magnitude. This shows that the 
kinematical dependence of rescattering functions is relevant in constructing the imaginary 
parts of the amplitudes. Nevertheless, such a large imaginary contribution to the term linear 
in Y does not help in generating the large CP-violating interference between the two 1=1 
amplitudes suggested in Ref. pl| . Indeed, of the two F-dependent terms in Eq. 



3), the 

one proportional to a c has the same weak phase as the constant term and consequently, as 
already noticed in Ref. PJlQl, the CP- violating interference between the amplitudes to the 
two 1 = 1 states must be proportional to the small difference (ao — A))- 

In principle, the rescattering phases should be included in the analysis of CP conserving 
Dalitz plot parameters. Their contribution could affect the determination of the linear and 
the quadratic slopes. However, since in this case the imaginary parts appear quadratically, 
their effect is of order p 8 in ChPT and thus for completeness also the other contributions of 
the same order should be included. As a curiosity, we estimate the contribution of 5smAQ 0Q 
to the quadratic slope h, in the Dalitz plot of K L — ► 3n°, defined by 



\A L I 2 
l^oool 



(JfteA 



000 ) i 



^smA L 



ooo > 



oc 1 + h(Y 2 + X 2 /3) + 



Using Eqs. (fL3D and (|28|) , we find: 

3mAo 00 = 3a n a + 3(b n (3' + a n ai)(Y 2 + X 2 /3), 
which gives the contribution to 



/i (Im) = 2a ( oti + —(3'A ~ +1.4 x 10- 3 . 



This number turns out to be of the same order of the experimental value of h 

h = -(3.3 ± 1.1) x 10~ 3 , 
but is substantially smaller than the p 6 contribution theoretically estimated in Ref. 



(31) 



(32) 



(33) 



(34) 



V. MEASUREMENTS OF THE RESCATTERING MATRIX IN 
INTERFEROMETRY MACHINES 

As pointed out in sect. |, measurements of Ki — K$ interference as a function of time 
should represent a convenient means to determine the (3tt) rescattering matrix elements, 



1 This numerical result is obtained with the value of b n /a n resulting from the fit of Ref. 
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because this observable depends linearly on 9fm[(A£_ )M£_ ]. To this purpose, "interfer- 
ometry machines" such as DA$NE and LEAR should have the advantage that interference 
naturally occurs in vacuum there. It is possible to measure Kl — K$ interference terms also 
in fixed-target experiments where statistics can be higher, however in this case an accurate 
knowledge of the regeneration amplitude is required. 



Recent LEAR data [18] give a preliminary indication of the term proportional to 
cos(Ami) in Eq. (1) and suggest the possibility of measuring in the near future also the 
sin(Amt) component. Consequently, it is worthwhile to improve the order of magnitude 
estimates of Ref. JT2] and, using the ChPT results of sects. [TIT], [TV], to derive predictions for 
machines like LEAR and DA$NE based on that definite theoretical model. 

Choosing \K°) = CP\K°), the CP-even and CP-odd eigenstates are \K Xy2 ) = (\K°) ± 
\~K°)) /V%, and, with the Wu-Yang phase convention, the mass eigenstates are (assuming 
CPT invar iance) 



\K s , L )= P \K°)± q \W)= ^ ) + £lK ^ 



1 + \e\ 



(35) 



The proper time evolution of initial K° or K° states is 
\K°(t)) 



V2(l 



\K s )exp 



-r s t 



- im s t + \K L ) exp 



-Trt 



— im^t 



\K°{t)) 



i + kl 



V2(l 



\K S ) exp ( — — ^ im s t) — \K L ) exp ( — — im L t 



(36) 



At LEAR, tagged K° and K° are produced, and the simplest means to observe interfer- 
ence is represented by the asymmetry 



/d$ f(X,Y) 


)A(K° -> TT + TT-TT )! 2 " \MK° -> TT + TT-TT )] 2 " 




'\A(K° -> TT + TT-TT )! 2 + \A(K° -> TT + TT-TT )! 2 " 





(37) 



where f(X, Y) is an odd- A function chosen in order to disentangle the different kinematical 
dependences. Up to first order in e, the decay amplitude squared as a function of time is 
given by 

\A (K\W) - Tr+Tr-Tr ) | 2 ~ ~ (1 T 2^ee) { exp (-T s t)\A s \ 2 + exp (-T L t)\A L \ 2 

±2 exp (-Ft) [3?e (A L A* S ) cos (Ami) + 3m (A L A* S ) sin (Amt)} }, (38) 

where Am = mi — ms, T = (Tl + Ts , )/2 and As,l = ^.+'-q- Then Eq. ( p7|) can be rewritten 



as 



2e- rt J d§ f(X, Y) lkeA L $leA s 



where 



rf$ [e- r st\A s \ 2 + e- r ^\A L \ 2 ] 



cos (Amt) + S f sin (Amt) + 0{8) 



(39) 
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J d<$> f(X, Y) [<3mA L MeA s - %mA s MeA L ] 
/d$ f(X,Y) MeA L Me A s 



(40) 



At the planned DA$NE machine, a Ks — Kl coherent state will be produced and the 



interference term of Eq. (ffl) can be studied by looking at the final state (/+ 



7l ± l>, 7T + 7T 



-7T°) 



Following Ref. [*2"""1 , we define for a generic decay K s ,lKl,s ~ > fiitijfiit?) an intensity 



/(/i,/ 2 ;t) 



1 /"°° 

2 7|t| 



rfT|</i(*i)/ 2 (* 2 )|*>| s 



(41) 



where t = ti—t 2 and T = £i + £ 2 - Choosing /i = Z+tt 11 // and / 2 = 7r+7r t 1 " , we can define an 
asymmetry similar to A~t~ (t), namely 



Rf(t) 

Indeed, for t > we have 



/ d$ f(X } Y) HF^u, tt+tt-tt ; t) 
Jd<$> /(Z+tt+z/, tt+tt-tt ; t) 



(42) 



/(Z+tt^, tt+tt-tt ; t > 0) = { exp (-T s t)\A L \ 2 + exp (-T L t)\A s \ 2 

±2 exp (— Tt) [Me {AiA* s ) cos (Ami) — Qtb (A L A*)sin (Ami)] }, 

and therefore 



> 0) = ± 



2e~ rt J rf$ /(X, Y) MeA L Me A. 



Jd$ [e~ 



rv.ti 



[cos (Ami) — 5/ sin (Amt)] 



(43) 



(44) 



where 5/ is the same as defined in Eq. ( "4"0"D and terms of order 5 2 have been neglected. For 
t < the analogue of Eq. ( 4*4]) is 



Rf(t < 0) 



2e- r l*l / rf$ /(X, F) KeA L MeA s 
~ fd<5> [e^sl^Asl 2 + e- r M\A L \ 2 } 



[cos(Am\t\) + 5 f sm(Am\t\)]. (45) 



quickly becomes much 



However, due to the exchange Ti «-> Ts, the denominator in Eq. ( 
larger than in Eq. ("44]), and suppresses the interference effect. 

Considering for 5/ the first non- vanishing order in ChPT, which is 0(p 6 ) in the numerator 
and 0(p 4 ) in the denominator, we obtain 
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J d§ f(X, Y) a n {a - a' - S)X - b n (P - - S)XY 



f 



(46) 



fd$f(X,Y) [a n X-b n XY] 

If we use as weight function f(X, Y) = sgn(Jf), we obtain numerically the following result: 

5 X = 0.18 ±0.01. (47) 

Essentially, this turns out to be: 5x = «o — ^o, and is practically independent of the 
theoretical uncertainties on the small ratio b n /a n . For this reason the result ("47|) is in good 
agreement with the prediction of Ref. [T~2[ . 
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On the other hand, choosing /(X, Y) = sgn(FX) we obtain numerically: 



S XY = 0.30 ± 0.05. (48) 



This result is about a factor four larger than obtained in fL2"f, Sxy — 0.07. Indeed, by 



expanding the rescattering functions, in the present calculation we have 

(/?o-^o) + T a K-25 1 ) 

Sxy * • (49) 

/ d$ |X|sgn(Y) 
hi. J. 



d$ \XY\ 



Eq. ( |49"D shows that also in the case of Sxy the Y-dependent terms give a sizable contribution, 
since they are multiplied by the large factor (a n /b n ). The error in Eq. ([|8|) accounts for 
the theoretical uncertainty on (a n /b n ), for which either the experimental value or the 0(p 2 ) 
ChPT prediction can be used.Q 



VI. CONCLUDING REMARKS 

In the previous sections we have introduced a general formalism to consistently account 
for final states interactions in K — > 3ir amplitudes, and have used leading order Chiral Per- 
turbation Theory to evaluate the rescattering matrix. We have considered some potentially 
observable effects of rescattering on Dalitz plot variables. The results indicate that the off- 
diagonal elements of the rescattering matrix in the 1 = 1 sector induce sizable imaginary 
parts in the X and Y dependent amplitudes. However, these large imaginary parts are not 
easily detected from Dalitz plot analyses and cancel in direct CP-violating asymmetries. 

Planned experiments at "interferometry machines" can have direct access to the rescat- 
tering matrix elements via appropriately defined time-dependent asymmetries, which we 
have estimated in leading order ChPT. As examples of the typical effects expected in this 
framework, Fig. 2 shows the asymmetry A x ~°(t) of Eq. (^) relevant to LEAR. The solid 
line represents the asymmetry with no rescattering (Sx = 0), the dashed line corresponds to 
the leading ChPT estimate of Eq. (f4~7|) , and finally the dotted line would result by doubling 
the value of Sx- To obtain Fig. 2, for the real parts of the amplitudes A+_ and A+_ we 
have used the expansion (|ITD with the values of the parameters obtained in the fit of 0. As 
one can see, the curves in this figure have similar shapes, but possibly could be distinguished 
in high precision experiments. 

In Fig. 3 we report the asymmetry R x (t) of Eq. ( [4"4] ) relevant to DA$NE, and the three 
curves refer to the same cases considered in Fig. 2. Here we note that rescattering affects 
the shape of the curves more significantly, especially for t > where the asymmetry can 
become quite large. However, this occurs for the values of t where the number of events 



2 The unknown quadratic term in ^4+_o can affect to some extent the numerical result for 5xy H > 
and its effect can be roughly taken into account by doubling the error in Eq. 
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becomes smaller. As an indication, the total expected number of events at t > 0, with the 
planned DA$NE luminosity 5 x 10 32 cm~ 2 sec^ 1 , is of the order of 10 3 /year. 

In conclusion, the analysis presented here shows the interest of experimental efforts to 
accurately measure the kind of asymmetries discussed here. The ultimate goal would be the 
determination of the K — > 3ir rescattering matrix elements testing ChPT in the strong sector, 
but in any case even a reasonable upper bound would represent an important information in 
this regard. Furthermore, the direct measurement of the CP conserving K$ — > Sir amplitude 
is by itself an important achievement, extremely useful in order to test chiral symmetry in 
non-leptonic weak interactions. 
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FIGURES 



FIG. 1. Loop diagrams relevant to K — > 3tt rescattering. The symbols • and o indicate the 
weak and the strong vertices, respectively. 

FIG. 2. The asymmetry A x ~° of Eq. (39) vs. t. The full, dashed and dotted lines correspond 
to 5x = 0, 0.2 and 0.4, respectively. 



FIG. 3. The asymmetry R x of Eq. ( |44"D for positive and negative t. The full, dashed and 
dotted lines correspond to 5x = 0, 0.2 and 0.4, respectively. 
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